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ABSTRACT 


One cof thevgreatest problems faced «by designers of 
simulation languages is producing a language whereby a 
user can construct a model which adequately represents 
the real system being simulated. This thesis will review 
the problem’ and illustrate the attempts made to solve it 
by reference to existing simulation languages. 

Two new simulation programs will then be described 
and si llustratedyeeathe first, running in the conventional 
batch stream, presents a simplified approach to the group- 
ing of model variables. The second represents a new 
approacheto the construct voneotesumulatiion models as it 


operates within a time-sharing environment. 
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CHAPTER I 


INTRODUCTION 


Lips Background 


One consequence of the development of the modern 
computer is that simulation has become an important 
technique for solving problems in many areas. These 
areas range from the study of the smallest molecules 
to the.most complex interactions of man and space. 

Although simulation has proved to be an important 
and valuable tool for use with a computer, it ae ae a 
product of twentieth century technology. For example, 
Lewis (1967) presents examples of simulation experiments 
conducted in each of the 17th, 18th, and 19th centuries. 
However, application of the technique to study the 
dynamic behaviour of large and complex physical systems 
did not occur until the development of mathematical 
techniques such as Monte Carlo Methods (Hammersly and 
Handscomb (1964)). Although use was made of these 
methods before the advent of computers (Farmer and Coll- 
cut (1967)), the feasibility of their use was greatly 
increased by the availability of computers. 

As indicated by Wallace (1967) simulation, as we 
Siglrecetinesit,.nad its origins In the study of. systems, 


a term which may be loosely defined as a collection of 
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parts capable of interacting in a manner that satisfies 
a specific set of requirements. The system may be real 
or hypothetical, and may be composed of components or 
Subsystems which may be permanent or temporary members 
of- the system. At any given instant the condition of 
a system is determined by the state of its subsystems. 
Thus, the state of a system at any two instants is the 
same.if and only if the state of all components of the 
system is the same at the two instants considered. A 
State history of a system may be defined. in terms of a 
given time interval and is the set-of conditions of the 
System at each of a set of successive time intervals. A 
partial state history of the system is obtained by col- 
lecting the conditions of the system at selected instants 
or by choosing the conditions of selected subsystems. A 
model is a representation of the system under study. It 
may be a physical model, mathematical model, or symbolic 
model. The state history obtained for the model is 
regarded as a state history of the system. 

We can then define the art of simulation as the con- 
struction of a model of a system and the collection of a 


state history of the model. 
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dae Early Developments of Simulation Models 


One of the earliest developments in the subject of 
simulation was the separation of simulation models into 
two distinct groups, continuous change and discrete 


change models. 


1.2.1 Continuous Change Models For many years the 


differential equations derived in the physical sciences 

and engineering have been solved by use of analog computers. 
Thus, systems which may be modelled by suitable differen- 
tial equations may be simulated by analog computation. 
However, problems associated with analog equipment include 
scaling of variables and non-linearity of equations. This 
has led to the development of digital analog simulators, 
Le programs) written-for, digital icomputers; but, which 
provide the basic operations of analog computers such as 
SsummMers-.. Lnjerrators:,.multip liens, .etc. Excelent) com— 
parisons of various continuous change languages for digital 
computers may be found in Brennan and Linebarger (1964, 


1965), and Clancy and Fineberg (1965). 


1.2.2 Discrete Change Models The present investiga- 
tion will deal with the second type of simulation model, 
the discrete change model. Teichroew and Lubin (1966) 
suggest that they may be characterized by the following 


properties: 
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1. The system is composed of subsystems of com- 
ponents, each of which has a certain prescribed 
mincviton tomperform. 

2. Certain items flow through the system from one 
subsystem to the next. Each item remains at a 
specific subsystem until its function is com- 
pleted; it ais) then transferred: tothe next 
Subsystem. 

3. Each subsystem has a capacity which cannot be 
exceeded. Thus, items may have to form a queue 


before being accepted by a subsystem. 


ihe first programs that utilized computers for simu- 
lation studies were written. to simulate one specific 
system in the language currently implemented on the com- 
puter. This language was often an assembler language. 
With the advent of compiler languages, such as FORTRAN, 
the programming was somewhat simplified. However, since 
the development of simulation models often involved many 
changes and consequent reprogramming, it soon became 
apparent that the special purpose models were not suf- 
PLletenvily flexible. 

This led to the development of general purpose simu- 
lation languages. At first these languages were very 


restrictive in their field of applicability. For example, 
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some: were ideveltoped particularly for job shop -simulations 
(IBM (1960)), inventory management simulations, and mili- 
tary operations. Then in 1959 one of the first of the 
true general simulation languages, MONTECODE, was intro- 
duced in England (Kelley and Buxton (1962)). Since then 
many general simulation languages have appeared and, as 
would be expected, each one has offered some different 
approach to the problem. 

As indicated by Kiviat (1967), with the appearance 
of general purpose languages simulation theory was soon 
divided into two schools, the flow chart approach of 
GPSS (Gordon (1962)) and the statement language approach 
of SIMSCRIPT (Markowitz, Hausner and Karr (1963)). 

The GPSS user specifies the simulation model by 
constructing a block diagram of the model. Each block 
is: selected from an.assortment provided by GPSS, and 
only these blocks may be used. The user does not specify 
exactly what happens in each block, but-provides certain 
parameters which act to specialize the.action of each 
block. 

The: models constructed from SIMSCRIPT, while often 
based on a flow diagram, are built up using program state- 
ments as in other higher level languages. Although the 
flow chart approach is easier to learn it is not nearly 


as flexible as the statement approach. For this reason 
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most simulation languages use the latter approach. 


Po aL Urposesofy the Study 


In the few years since the introduction of general 
purpose simulation languages, a start has been made to 
develop a. theory of simulation. The present thesis will 
review the theory which has developed and will illustrate 
the important points with reference to existing languages. 

Two new. simulation programs will then be discussed. 
The soLrpst tsoewriccenein FORTRAN ang@eis i desiened: to: run in 
normal bateh mode. = The second iAs-written in APL, an 
automated version of a notation which was first described 
in Iverson's "A Programming Language" (1962). It is im- 
plemented in a time sharing environment and represents a 


new approach to the development of simulation programs. 
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CHAPTER II 


THEORY OF DISCRETE CHANGE SIMULATION 


ea I Requirements Of 4a Simulation Language 


A simulation language can be distinguished from 
other computer languages by the special features or 
services it must. provide to a user who writes a computer 
program based on a model of a system. Numerous authors 
have presented their estimates of what form these dis- 
tinguishing services should take (Krasnow and Merikallio 
(1964), Farmer and Collcut, (1967), Krasnow (1967), 
Teichroew and Lubin (1966), Jones (1967a), and Kiviat 
(1967)). The ones most frequently mentioned are listed 


below. 


1. Since. the model.is represented in the, computer 
as a set of numbers which describe the state of 
the model at a particular instant, the language 
must provide a means whereby this data set can 
be structured into variables, data, queues, 
etc., each corresponding to some recognizable 
entity of the system. being simulated. For 
example, SIMSCRIPT allows a description of the 
state of a model in terms of entities, attributes 


Of ett Ges] and ereups.of entities called sets. 
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Krasnow and Merikallio (1964) refer to these 
as status descriptors. 

The simulation language must contain certain 
instructions with which the programmer can 
change the state of the model. These include 
instructions for moving members of a set, 
arrenmetic Instructions, and instructions) to 
change values of. the data base. 

The sequence in which state changes are to 
occur is often highly complex and unpredict- 
able. pThus,. each» language, should contain, a 
built-in timing and scheduling mechanism. 
Facilities to provide statistical and, other 
performance data, and to produce meaningful 
output. 

An important requirement of any computer lan- 
guage is the debugging facilities. These are 
especially important in simulation programs 
since use of random number generators and 
probability distributions make reproducible 
runs hicult. to obtain, 

Facilities for generating random numbers and 
for sampling probability distributions. GPSS 
contains excellent facilities for user and/or 


system supplied functions. 
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It must be noted that the above six features can- 
not be regarded as exhaustive or necessary in every 
language. However, the manner in which the language 
implements: these features is important and this will 


be. discussed-aintthe  Tollowing’® sections 


2.2 World View 

The world view of a simulation language is the con- 
ceptual foundation upon which the system may be modelled 
(Krasnow and.Merikallio (1964)). It is the way in which 
the language arbitrarily divides the resources and con- 
cepts. Of thet real@ worid: toe provide the’ services’ listed 
in- Sectiono2.1°°°To the user-it represents the way in 
which he must think about his- problem in order to use 
the language to develop a simulation model. 

In addition to providing the services previously 
mentioned, the language must resolve certain aspects of 
the real world and take them into account in its world 
view. Krasnow (1967) lists the following four character- 
istics of the real world the simulation language must pay 


particular attention to. 


1, Continuity’ The» real world is made up of con- 
tinuous actions while the digital computer is 


capable of producing only discrete changes. 
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2. Uniqueness It would be impossible to write a 
program describing every unique action of the 
real world. It is necessary to isolate basic 
repeatable patterns and use these as repeatable 
units in-the=simulation "program, 

3. *Interdépendéncé* Itristimpossiblestotaccount for 
the mutual interactions of all activities or 
events in the real world. The language designers 
must then determine how much interaction of the 
real world will be modelled. 

4, Simultaneity At any instant in the real world 
there occur infinitely many events. The digital 
computer, however, can perform one. or, at most, 


a few operations at one time. 


The remainder of.this section is devoted to a dis-= 
cussion of some of the basic concepts which are used in 
describing the world view of a simulation language. It 
is based on papers by Krasnow (1967) and Kiviat (1967) 
which contain comprehensive discussions of the concepts. 

The writing of computer. programs.as subroutines or 
procedures has proved to-be an efficient and meaningful 
method. Although actions and events are unique in the 
real world it has been found essential to provide the 


same facility in simulation languages. This facility is 
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usedebyhrecognizing certaintpatterns cof «behaviourain 
the real world and differentiating each occurrence by 
specifying parameters. 

Torillustrate the use of patterns of behaviour 
consider automobiles entering a service station for 
gasoline. The pattern of car entering, possibly wait- 
ing in line, getting served, and leaving is repeated 
for léveny car. 

tere patitennican theewritteneassacsubreutine: or 
set of subroutines with parameters specifying gasoline 
type and quantity. These basic modules are referred to 
as Mdesceiptive sunits het athe elanguage:, landarepresent: the 
repeatable patterns of some activity in: the system being 
mode bled:.esgin-activity. refers sto etheideseriptioneof the 
state-«changes: whichoccur withinvardescriptive unit: 

We fare enow facediewith rthesprobiem cof ‘resolving’ the 
concep Girof Jaicdescriptive! unishitwilth atheefact ofsa con-— 
tinuously changing and interacting real-world. Two 
related considerations may be mentioned.» First, in:the 
real world novaction eceurs. without the passage -of time. 
However, a computer is capable of producing one, state 
change at a time, and it is accomplished without the 
passage of simulated etime:.i | Second:.,rar state schange scaused 
by,.one, activity may influence. the: action »of another activ- 


ity, to.perfonm another, state change. 
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To check for all impldeattions after a state change 
would add considerably to the programming burden of the 
user and implementer, sincevxa descriptive unit: could 
contain only one state change. The problem is resolved 
by making a number of state changes within a descriptive 
unit, all at the same simulation time, until a plausible 
interaction point is reached. 

iivalleche stave changes withingan activity eccur 
at the same simulation time, and terminate with-an inter- 
SclleOnepoint an execution er instance of such an activity 
is an event. In 1964 with the announcement of SOL (Knuth 
and McNeley (1964)), the concept Oeftaneinstance ef-an 
activity was extended to include activities in which the 
Suave cnanges~oceur over a period of=simulated time., An 
Pisvencemol suchvangactivity sea sonocess,.) Then 1f tan 
activity has only a single interaction point, event and 
process are identical, and event is a special case of a 
process. Note. that in the above - discussion activity 
refers to the’ state changes which can occur while event 


and: process refer to the execution of the activity. 


2.3 Machine-Based vs. Material-Based World View 
The larger differences in existing simulation lan- 
guages can be traced to the type of descriptive unit used 


and how they are used to develop a world view. Two separate 
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approaches can be recognized, the machine-based and 
material=-based approach. Each will be discussed with 


examples of existing languages illustrating each approach. 


2.3.1 Machine-Based Approach In this approach the 


language allows the programmer to segment the program into 
activities (usually subroutines) which contain instructions 
to make state changes. At every simulation time an attempt 
is made to execute each of the activities in turn. The 
execution of the state changes in an activity is controlled 
by a set of test statements located at the beginning of 
each- activity. Jf all-conditions are met the rest of the 
activity is executed, otherwise the next activity is im- 
mediately attempted. After no further activities can be 
executed at a simulation time, the clock is updated and 

the procedure is repeated. Note that this type of organ- 
ization requires no separate scheduling algorithm aside 
from a control program which branches to each activity in 
turn. 

Updating the clock can be handled in either of two 
ways; using a fixed time increment or advancing the clock 
to the time the next event is to occur. The latter approach 
ae ayreok provide a more realistic simulation (Evans, 
Wallace and Sutherland (1967)) but will require a schedul- 
ing algorithm which scans a list of event notices to find 


the next event time. 
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2.3.1.1 Control and Simulation Language CSL 
(Buxton and Laski (1962)) is an example of a machine-based 


language which advances simulated time by scheduling 
events. The language provides the programmer with special 
status descriptors called T-cells, which are a special 
type of integer variable. If, for example, the programmer 
required a variable SHCHANGE to account for a shift change 
in a simulation, he would define it by specifying TIME 
SHCHANGE, and in an activity it would be referenced by 
T.SHCHANGE. 

The function of: the T-cells in updating the clock 
Can be déseribed as follows. After all activities have 
been executed (or attempted) at a given simulated time, a 
scan is performed on all.positive valued T-cells to find 
Chesmitnamun. wALISl—cetis are reduced by the minimum 
value and the clock is advanced by this value. Thus, at 
teast one T-cell will become zero and, during the. next 
execuvioneoleuhe activities, if an activity only centains 
a, test for this T-cell = zero, the activity will be com- 


pletely “executed at this updated clock time. 


2.3.2 Material-Based Approach Two methods have 


been used to implement the material-based approach, the 
event orientation of SIMSCRIPT and the process orientation 
of SOL and SIMULA (Dahl and Nygaard (1966)). Although the 


event was previously described as a special case of the 
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process, examples of each will be described since the 
languages resulting from the two approaches are quite 
different. Regardless of the orientation used, the 
material-based approach differs from the machine-based 
in requiring a control algorithm which scans over a 
schedule of future events to determine the next state 
changes (note that in the machine-based language CSL a 
scan was made only to determine the next clock time, not 
which activity was to be executed next). Since each 
future event must have a time associated with it, this 


scan also determines the next clock time. 


2.3.2.1 SIMSCRIPT The purest event oriented 
language is SIMSCRIPT. The programmer constructs a set 
of descriptive units called EVENT ROUTINES (similar to 
FORTRAN subroutines), the execution of each constituting 
an event as previously defined. An event is scheduled 
DyetosulneeanecVetu nNoulce ,.and specifying ai time at 
Wit GmeLnereVenl isto occur.) “hiss is accomplished: by 
the statement CAUSE (event notice) AT (time). After the 
execution of an event routine control passes back to the 
tim~nesomecOntrolsroutine which scans the list of scheduled 
Eventenovlces slo Pind the next scheduled event. The clock 
Pemupdaved watias eEvenu novice is deleted from the list, and 
the event routine associated with this event notice is executed. 


The procedure continues until no event notices are scheduled. 
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2.3.2.2 Simulation-Oriented Language The 
introduction of SOL in 1964 provided simulation program- 


mers with a new approach in the construction of simulation 
programs. A SOL program consists of a set of descriptive 
units which in our notation can be classed as activities 
(the authors refer to them as processes). The activities 
describe the actions performed by certain entities of the 
real world, with the programmer having panieremeontna! 
over the number, definition, and detail of the activities. 

For example, to simulate traffic flow, a SOL program 
might consist of one activity describing the movement of 
automobiles, another describing the movement of trucks and 
buses, another describing the action of traffic lights, 
and another describing the movement of pedestrians. All 
activities are considered to be operating in "parallel" 
with all others, with interaction being achieved through 
changing of global variables which are accessible by 
every activity. 

AS atwGros , ~schedusine, and timing in) SOL are aided 
DyMuUScHOlepransacti ons which) can be used to» represent the 
inch nOimo Venporary.cnuilty.s through) the model, Since, the 
Lroanscacvichmiscmune OCcCurrence Of an activity, it can be 
called a process. Then, for example, a transaction moving 
through the automobile activity in the previous example 


would represent an automobile moving through the model. 
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There may exist any number of transactions in a 
many-to-one relationship with an activity, each executing 
different statements within the activity. Each trans- 
action has associated with it a time variable, a pointer 
to whére the transaction is located in the activity, and 


local variables. 


2.5-c.5 wolMUlation LAnguage The designers of 


SIMULA recognized the process concept of simulation lan- 
guages as similar to procedures of algorithmic languages 
such as ALGOL (Naur (1963)), and consequently extended an 
ALGOL 60 compiler to allow simulation programs to be 
written. One of the major changes necessary was to allow 
"narallel" execution of the simulation activities instead 
of sequential execution as in the usual implementation. 
As with SOL the model is written as a set of activi- 
ties; tne execution of an activity being a process. Once 
a process has been initiated it may be in one of four 
states; active, passive, suspended, or terminated. Only 
one process can be active at one time, corresponding to 
the movement of only one transaction in SOL. When a pro- 
cess reaches a point in an activity which causes its state 
to change from active to passive or suspended, the control 
program automatically issues an event notice for that 
process, associates a time for next event with it, and 


Dicceseiy Lieene sequence seu, a schedule of next: events. 
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It also associates a reactivation point which specifies 
the location within the activity at which the process is 


to begin executing when it next becomes active. 


2.4 Summary 


As is evident from the preceding discussion there 
exists considerable variety in the approaches used to 
nesolve the dynamic laspects of the real .world., Of<those 
developed so far the process approach of SOL and SIMULA 
appears to .be the most flexible, and in a recent paper 
by Blunden and Krasnow (1967) the authors propose the 
process concept as a basis for a new simulation language. 
However, as indicated by Teichroew, et al (1967), 'there 
is as yet no conclusive evidence that one simulation 
language is best for a variety of problems'. 

There are a number of papers which compare the 
languages mentioned in thischapter (and others) and the 
approach used by each in detail. See Tocher (1965), 
Krasnow (1967), Teichroew and Lubin (1966), and Krasnow 
and Merikallio (1964), all of which also contain excel- 
lent bibliographies. For a detailed description of a 
particular language, see the reference manuals available 


for that language. 
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CHAPTERS DEL 


A MACHINE-BASED SIMULATION PROGRAM 


Se lneroductilon 

This chapter will describe a simulation program 
which uses the machine-based world view as described 
in Section 2.3.1. The program is written in System 360 
FORTRAN IV source language and a complete listing has 
been included in Appendix I. 

It is supposed that the system to be simulated can 
be represented bysa flow dlagram consistdrigsof two basic 
components, blocks of entities and activities. A seg- 
ment of a typical flow diagram, showing the two compon- 
ents, is given in’ Kigure 1. 

The blocks. of entities are used to record the state 
Of variables required in the simulation model, at. various 
locations in the system being modelled. The position of 
the variables on the flow diagram corresponds directly 
to the position of the entities they represent in the 
Rea uesyc ce selOor example. sah i oure ss. , ther first entity 
in each block may represent the number of automobiles 
Dresemoe au that particular docataon in the model, while 
the second entity may be used to represent the state of 
trabpiculichossaterhe same locatlons. Changes in-the 


state of the model are recorded as changes in the values 
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Block 1 


Block 3 Block 4 


Activity 1l 


Activity 2 


Block 2 


FIGURE 1 


eal 


of- the state entities’ in some, or all, of:the entity 
blocks. The user must assign the function of each 
entity required in a model as it is being developed. 

As indicated in Figure 1 the blocks of entities 
are joined to activities. An activity contains, prim- 
arily, instructions to change the state of various 
entities. The programmer indicates which entities are 
to be changed by linking those entities with an activity 
which will cause the desired state changes. After the 
flow diagram has been constructed each activity is 
written as a FORTRAN subroutine subprogram. Details of 
their construction will be presented in’ a Tater’ section. 

simulation time is recorded in the variable TIME 
and is available to all activities. Time is advanced by 


means of a fixed time step supplied by the user. 


3-e Details of the Program 


7 enoGliees andevanclanleswmeAssindicated in the 
previous chapter any system to be simulated contains 
certain entities and/or variables which describe the 
state of the system at a particular time. In the present 
Procraneat a Veriablessand entities required to define: the 
system must be in one of the entity blocks. The entity 
blocks are numbered sequentially beginning at one, and 


the individual entities within a block are numbered in 
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the same way. The entities required are arranged as a 
two-dimensional array referenced by the programmer as 
STATE(I,J). The first dimension refers to the entity 
block number, the second to the entity within that 
block. The number of entity blocks and entities per 


block are defined by the user, as will be described later. 


Beewe ACUlViLties 


Bivie eeiek Configuration Parameters As described 


in Section 3.1 each activity may have access to only 
certain specified entity blocks. The user can write the 
activities in one of two ways to specify the entities and 
entity blocks available to an activity. In the first 
method the user specifies explicitly which entity is of 
interest. For example, if an activity required the fourth 
entity in entity block two, the user would reference it 

ty the activivy as SlATE(c,4). 

In the second method the user specifies the entity 
references as integer variables. The entity references 
are then known as configuration parameters which are 
specified at the time the activities are being assembled. 
The activities can then be written so that they can be 
used at different places within the model, with each 
application differentiated by the configuration para-= 


meters. 
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AE cote Activity Parameters The user has one 


additional way to specify the uniaue action OT saneacuLvLcy 
This is by use of activity parameters. Whereas the con- 
figuration parameters are used to indicate the relation- 
Snip between entity blocks and activities, the activity 
parameters are used to differentiate between the action 
Of sune  salewaclivities. Or example, an activity may be 
written which will move a multiple of an entity from one 
block to another. This subroutine may be used a number 
of times within the flow diagram. The difference in each 
easels the valuevof ithe multiplication factor.) Both 

the configuration and activity parameters are supplied 
during the assembling of the model and remain constant 


for that model. 


Sericl SEO LMaUsOleACo AEDES ln eCOnsSurUuCccing 
the activities it is necessary to impose a number of 
restrictions over and above those which result from using 
the FORTRAN language. They pertain primarily to the use 
and definition of some variables. Figure 2 contains an 
example of a typical activity. The following outlines 
those features which must be present in each activity to 
be run using the present program. 

1. The subroutine name must be of the form OPNX, 
where, in the present program, xX is an integer 


from one to sixteen. This allows for up to 
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SUBRUOUTING UPNI 

CUMMUN/BLAA/ STATES PPAR pKRPAR,NAMUP, NKPARgNPPAK yg NLINUs INUMO 
CUMMUN/BLKGE/NOLUKyNSTAT,NKP ap NP Re KKP yp KRPP,yLNUEA 
CUMMUN/ SLKLS NCUINY NUPS »>MAG,NUSE » MUN yNUUT, NUAT ’ Tlie 
DIMENSLUN NAMGP( 16), KiNUL 144), NUME L144) ,RPAR( LOU)» 
LNKPAR(LLi6) sNPPAR(1i0), STATE( 25440) sPPAK(i00; 

NKP=3 

NPP=0 

GU TU t1y2),1NDEX 

RETURN 

I = KPARCKKP) 

J KPAR(KKP + 1) 

K = KPARKCKKAP + 2: 

Pe SAT las SU ede oO). KE.TLURN 

SUALEC Laws <= chaleliscKhy * DD 

RETURN 

END 
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sixteen activities to be defined and used in a 
model. The name of each new subroutine used must 
Demerniered Roy ta tcalMs tevaktspecial sectontor ithe 
control program, to be described later. 

2. The COMMON and DIMENSION statements are required 
to allow the subroutine and control program 
access to the data base. 

Jal Themexp tour, statements must appear asvindicated 
in the example, immediately following the DIMEN- 
SION statement. The two assignment statements 
define the number of configuration parameters 
(NKP) and activity parameters (NPP), respectively, 
fous Tactivityes= Thespurposelof this) setivof 
statements will be described in a later section. 

4, The remainder of the activity may be considered 
vo berdividedeinto twotsections.  Lhewfirst 
contains test statements which must be satisfied 
before the state changes contained in the second 
section are performed. All the state changes 
indicated in the latter section must be capable 


of being executed at the same simulated time. 


3.3 Control Program and Command Words 


The program offers the user a number of functions 


which can be performed by use of command words. The 
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simulation program is run by specifying one of the command 
words which causes a branch to a certain section of the 
control program. The program performs the function called 
for after which it returns and requests another command 
word. 

inescontrolproeram also handles! all input, and output, 
unless it is data requested by an activity during a simula- 
tion run. The program produces output resembling that 
obtained at a remote terminal of a time-shared system. 
All input. statements are preceded by a line of output 
requesting the desired information. Immediately after 
Ghee pute Se reader cron Cardsmitetsmprintedsout, chus pro— 
ducing a well-defined record of all input to the program. 

The program can also be instructed to store a model 
currently in memory on magnetic tape. This stored model 
may then be recalled and its operation continued at some 
haters datex 

The function and operation of each command word is 


described in the remainder of this section. 


Soe eV heli romCOmnandmesmUused  vosepmepare 7a 
new magnetic tape by rewinding it and placing on it a 
single record of nineteen words consisting of -1, 0, O, 
Ge bilaniewwordsy lhe first zero indicates that no activities 
have been defined, the second indicates that no models have 


been stored on the tape, and the sixteen blank words contain 
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space for the names!of-up tocsixteen activities. This 
record forms a header record on the magnetic tape. As 
other records are added to the tape the header record is 
moved along so that it always remains the last record on 


the tape and provides a summary of the previous records. 


3.3.2 NEWACT This command is used to define new 
activities. Following input of the command word the 
program requests the names of the new operations, each 
containing four characters followed by a blank. The pro- 
gram then assigns each operation a unique operation number 
chosen sequentially from the first sixteen integers. The 
header record on the tape is updated, containing the names 


of all the defined activities. 


3.3.3 NEWMOD This command is used to set up a new 
simulation model. In response to the command the program 
requestsetwomparameters: » They first ,,\NSTAT wisethe number 
of state entities required in each block of entities. The 
second, NBLOK, is the number of blocks of state entities 
required in the model. 

TheSusecrathenwtistsathes@activities tosbeusedsin the 
model in the sequence in which they will be attempted by 
the scheduling program. To differentiate between multiple 
use of the same activity, each activity name must be fol- 


lowed by an integer. The program then establishes in 
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memory a vector KIND to indicate the order in which the 
defined activities are to be executed. The components 
of KIND are the activity numbers that were defined by 
the NEWACT command. Thus, the activities are identified 
by name for the user but by number for the computer. 

For example, assume the activities LINT, DIFF, and 
HYDV were entered as the first three activities by the 
NEWACT command. Then if the sequence of activities is 


input#as 
Dit PISHYDVI DiFK eS bINT weap hr 3s = 


WAESVectTOr KEND wini¥econsist of@25es, M2et 1 =e) 

After the’ activities and their sequence has been 
established, the program establishes two additional 
control vectors NKPAR(I), NPPAR(I), I=1,NTOT, where NTOT 
is the number of components in the vector KIND and repre- 
sents the number of activities used in the model. The 
Te component of NKPAR is the number of configuration 


parameters required by the 1°" activity, and of NPPAR the 


number of activity parameters required by the gee ACG Viney. 
To evaluate the parameters required in the above two 

vectors the program branches to each activity in turn as 

specified in the vector KIND. As shown in Figure 2 the 


first executable statements in each subroutine assign values 


to the common variables NKP and NPP. These represent, 
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respectively, the number of configuration parameters and 
activity parameters required in this activity. During 

this phase of using the program INDEX has been set to 

unity. Thus, the subroutine will be exited immediately 

and upon return to the control program the moh component 

of NKPAR and NPPAR will be set to NKP and NPP, respectively. 
This procedure is repeated until all components have been 
evaluated. 

After the above two vectors have been evaluated the 
program requests, as input on cards, NKPAR(I) configuration 
parameters for each activity I and NPPAR(I) activity para- 
meters for each activity I. The configuration parameters 
are stored in the vector KPAR and the activity parameters 
in’ the vector PPAR. At this point the model is completely 
specified and if all activities have been compiled the 


model may be run. 


3.3.4 STORE This command causes a check to be made 
as to whether a model is in memory available to be stored. 
If available the program will create a new record on tape 
immediately following the last stored model consisting 
of all information required to retrieve and begin running 
the model at a later time. The header record is updated 
indicating an additional model has been stored and placed 


at the end of the new model. 
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3.3.5 CALL This command causes a previously stored 
model to be retrieved from magnetic tape and placed in 
memory as the active model. Each stored model is referred 
to by a sequential integer, the first model referred to 
as, One. en tf yohes user specifies: retrieval of model J, 
the ane model stored on the tape will be obtained, unless 


J is greater than the number of stored models. If the 


latter occurs an error message will be output. 


3.3.6 RUN This command allows a model, which has 
been assembled and for which all subroutines are compiled, 
to be run. The program requests the time increment and 
lengLh ofeasimulation run. as anpuu data. lt athen, branches 
to each subroutine in the sequence indicated in the vector 
KIND... After each cycle through KIND the simulation,.clock 
is advanced by the time increment. supplied by the user. 
The procedure is repeated until time has been advanced 
the specified amount. 

After completing each cycle through the activities, 
and before the clock is updated, the control program 
branches to another subroutine, OUTPUT, which the user is 
also required cto supply. It is” through this™subroutine 
that the user obtains the required state history of the 
model currently being run. Since it is referenced at 
each simulated time it is possible to obtain a complete 


Stave ttisuory of the model. ~ This*would be “quite. valuable 
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in debugging a model. However, it can be used quite 
selectively as shown by the example in the next section, 
where the value of a few entities are obtained at the 
ence Olwasrun.s APcopy of the output subroutine used in 


the example can be found in Appendix I. 


Seaei) NIT When Che control program encounters 
this command word it immediately branches to a subroutine 
with the same name, where the user is required to supply 
a yes or*nowanswer Co CWwol questwvons im turn. The first 
Pequivesecs yes uC 1nitialiZzevallestare.entivies #to Zero 
and the second a yes. to initialize only some state entities, 
as specitied by the- user. "For the Latter the user is re- 
quired to supply a card for each entity to be assigned a 
non-zero value. Each card contains the block number, 
entity number within the block,! and the value of the 


entity, with a zero block number terminating the cards. 


3.3.8 END This command signifies the end of data 


for the simulation program, terminating its execution. 


3.4 An Example - Simulation of a Simple Three Machine 
Assembly System 


oe eco Cr iInUL Oliv Om RGOMleMIeure 3apnesents a 
flow diagram of .a system that will be used to illustrate 


the features of the simulation program presented in this 
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20 minutes/cycle Machine B 


1 component per 1 component per 
machine A cycle machine B cycle 


2 components/machine C cycle 


Machine C 15 minutes/cycle 
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chapter. It can be regarded as a sub-system of a larger 
production system which depends on the output of final 
product from machine] C. 

Basically the system consists of three facilities, 
machines A, B and.C. Machines A and B are started up at 
the same time and each produces one component every twenty 
minutes, for a total of two components every twenty minutes. 
There is, however, a 10% chance that any component may be 
rejected due to imperfections. Any rejects are discarded 
from the process. Components not rejected are placed in 
aesterace pDin. 

Machine C requires two components from the storage 
bin to produces a irinished productim taking 15 minutes to 
complete. If the storage bin contains fewer than two 


components machine C will remain idle. 


3.4.2 The Model Figure 4 presents the final entity 
block-activity flow diagram developed for this problem, 
Snceconsistseomesourscenuutcy blocks ~with a maximum-orf: Chree 
entities per block, and five activities. Only fiour sub-— 
routines are required since activities one and five use 
the same subroutine denoted in the flow diagram as SETF 
and compiléeédvas Subroutine OPNl. It requires, the specifica- 
tion of three configuration parameters to completely define 
tas stne=Othem three subroutines do not require any con- 


figuration oreactivity paramewers supplied prior to running 
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the model. A listing of the subroutines may be found in 
Appendix I 


The entities used in the model are defined as follows: 


STATE(1,1) 20 minute timer flag 
STATE(1,2) 20 minute timer 
STATE(2,1) component from machine A 
STATE(2,2) component from machine B 
STATE(2,3) total components rejected 


STATE(3,1) total components currently in storage 
bin 


STATE (342)* 915 “minute simer flag 
STATHC3,3) 15 minute timer 
STATE(4,1) total finished products from machine C 


STATE(4,2) total time machine C in operation. 


The following 1s ardescription of each activity, in 


the model and of the state changes caused by each. 


ie och eerie. tnreeecontigunationtparamevers Ly J’, 
Kerequested by subroutine SETFVare 1,71, 2 
Fespectively ia Shih re wlhis activity simply 
increases the contents of STATE(1,2) by five 
whenever STATE(1,1) = 1. Since the flag is 
Hever meohancean bails will occur on each cycle. 

2. RANDI This activity acts in conjunction with 


SETF1 Co provide the model with a 20 minute 
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Viner = eaWoienever s STATH Gl <2) < 20.n0 saction,is 
taken, otherwise the following occurs., STATE(1,2) 
is reset to zero and two random numbers, say RN1 
and RN2, are generated. If RNi < .9 set 
SHATHC2 a Je= 1sotherwisescTATE( 2.1) .= 0 .1=1,2. 
The action of this activity represents the pro- 
duction of two components every 20 minutes. The 
random numbers are used to determine if a compon- 
Cte comomre |e ChuOr NOL mo en ted Dy placanela 
zero or,one, respectively, in STATE(2,1), i=1,2. 
For every rejected component increment ‘STATE(2,3) 
by one. 

STORL This activity adds the completed and 
acceptable components, in STATE(2,1), i=1,2, 

to those already.in the storage bin, represented 
Acme lA er lame olATE erie) Lal eare. immediately 
meset. to zero. Note thatethis, activity can, be 
executed every cycle since nothing will be added 
to the storage bin unless components are produced 
by the previous activity RAND1. 

MACHT This activity represents the state changes 
asseecwvaved with machine C:and can be divided into 
two sections. 

(a) A check is made to determine if machine C 
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not busy, and go to part (b). Otherwise, 
enecks OTATE(3 8 )-etherl> minute timer, to 
determine if the 15 minute processing time 
nasmetapsede Pre STATE @3 53) *<er> no *further 
action occurs. Otherwise set STATE(3,2) = 
STATE(3,3) = 0 and increment the count of 
completed products by one. 
(by"*Check"=1fethe= storage sbin contains at *léast 
two components to utilize machine C. If 
STATE(3,1) < 2 machine C cannot be used. 
Otherwise set machine C to being used by 
STATE(3,2) = 1 and decrease the amount of 
components in the storage bin by two. 

Deo oie meme emu nvec requimed cCentioiraylen=epara— 
Meters ware=s, 2,730 sthewactivary will then 
inerease “STATE (3,3) "by five 1f STATHEC3 2) ="2; 
Thus the fifteen minute timer is used only if 
machine C is busy, and indicates the aatetne of 
processing time used on the current set of 


components. 


Selesenesultcs. Appendix 2 contains a listing of the 
output obtained when the example of this section was 
Simulaved. |OuLpUL Originating In the control pregram is 
preceded by C) while lines of output corresponding to user 


supplied data are preceded by U). Any output not prefixed 
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in this manner originates in the subroutines supplied by 
BHeRUSeTRs 

Since no activities or models have previously been 
defined the user SRI NEWTAPE to initialize the 
counters of these entities. The next command word 
specifies that new activities are to be defined. The 
program then assigns each a subroutine name by which the 
activities must be compiled. Then, with the activities 
compiled a new model is set up, containing four entity 
blocks with three entities per block. The next line of 
user data convains the names of the activities used in 
whe model in the sequence they are to be executed. Each 
name is followed by an integer indicating the number of 
tames the subrouvine hes been used up to this point. 
Next the configuration and activity parameters are input 
and ea nocrminlualiziie une enuLuves , the model is ™ready. 
Camvpe ful. 

After supplying the time increment and length of 
SLi ovLOWeG lhe uniesmOodel is ECXeCuvedmending wavhea 
Louie OmeuiesOuLpUL Of MachinesG, and the total time 
machine C was running. As indicated in the listing of 
the results in Appendix I, in an eight-hour run machine C 
produced 19 finished products and was in use 290 out of 


480 minutes. 
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CHAPTER IV 


OSS - AN ON-LINE SIMULATION SYSTEM 


Polen roductien 

ie beeneasonsetor an, On-Line System As indicated 
by Dimsdale and Markowitz (1964) and Jones (1967a) the 
normalecourse followed in simulation is: (1) writing of 
programs to model the system, (2) testing of programs, 
and (3) generation of output from simulation runs. The 
first two steps, however, are usually repeated many times 
in order to obtain an” accurate model of. the system, 

Conventional methods of writing simulation programs, 
using languages such as described in the previous chapter 
GEYGPSSPVSOIMSCRERPT, © SEMULAT eteoy require *the programmer 
ter SUbMIity his programs tO be run aim the off=line batch 
processing Stream. This method: usually implies a large 
number of-runs to obtain a valid and debugged model, 
resuLvuine wi nthe loss Of considerable vime waiving, for 
The=preegram vo be returned from the batch stream. 

A recent development in modern computing technology 
is the technique of time-sharing. Special purpose hard- 
ware and software systems allow a number of users at 
remote terminals to simultaneously have access to a large 
computer. 


On-line computing facilities provide means of 
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reducing the time required to build and test a simulation 
model since the building and testing may be done more or 
less simultaneously. One on-line system for simulation, 
OPS-3 (Greenberger, et al (1965)), has already been 
implemented and a new language and system for on-line use 
has recently been specified (Jones (1967a, 1967b)). 

This chapter will describe a new on-line simulation 
system written in APL and used on the IBM 360/67 at the 


Unversity sole Alberta. 


WoL ate she anuresieo 8 rantOn-line System As discussed 


in Section 2.1 all simulation languages must provide 
certain features if accurate models are to be produced. 
In an on-line environment a number of additional features 


may be considered desirable. Some are listed below. 


1. The system should allow the programmer to build 
CHomvecSumSeCUL ONS On thewmocel, peforesallythe 
PLeces have: been put togetner. 

2. An interpretive and compiled mode of execution 
should be available; interpretive for building 
and testing the model to avoid recompilations 
and compiled for running completed programs. 

3. Comprehensive on-line debugging and tracing 
facwi icles, should be vine luded. 


ie There should be facilities to save a model at 
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any =stave®or* construction -or=running "and recall 


it at any later time. 


eS. Deseription of OSS 


This section will present a detailed description of 
the OSS system. See Appendix II for a complete listing 


ofthe APL programs. 


4.2.1 General Organization and World View The world 
view of OSS encompasses the material-based or flow oriented 
models it ts assumed in this approach that the dynamic 
behaviour of a system is caused by something moving through 
the system which acts in some manner to cause changes in 
the state of the system. 

The state of a system modelled in OSS is indicated 
by two types of entities, permanent and temporary. Per- 
manent -entities"are present during the entire simulation 
whereas temporary entities are created and destroyed 
aderine?a simulation=run? = In°-O0Seethe*latter “are. catted 
Gransactions, and are used to represent those items which 
can be thought of as flowing through the system. 

AV@simulation ®prerram invOSS?consists “ofa set of 
aetivities@or “runetions in APL notation *(throughout “this 
discussion the terms will be used interchangeably). Each 
activity contains instructions for state changes and, 


Since the execution of an APL function cannot be suspended, 
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to continue at a later time, the state changes in an 
activity must all occur at the same simulated time. Thus, 
each execution of a function constitutes an event as in 
SEEMS CRIERT, 

The transactions flow through the model by causing 
une execut.onsot athesactivities, ateparticular apolints .1n 
time. OSS automatically maintains lists of when the 
events are to occur, and a special variable, referred 
HO, as. aclock iorssimulation time. peacords .the instant of 
real time that has been reached in the simulated model. 

AL) yctocketimes in OSS are -integens,«the unit of, 
real time corresponding to a unit of simulated time being 


defined by the user and used throughout a given model. 


4.2.2 Transactions Although the use of transactions 
is most apparent when simulating a system in which physical 
objects are flowing, in a more abstract view they may be 
interpreted only as vehicles for scheduling activities 
(Krasnow (1967)). Either view may be used when planning 
a model. 

hech transaceion in, OSS consists of a set: of integers 
divided aneo two classes; a group of six system variables 
used by the scheduling programs, and a variable number of 
user defined parameters. 

The system variables consist of a transaction priority 


number, a transaction number, next event indicators, a 
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time of next event indicator, and a transaction status 
indicator. The only one directly accessible to the user 
is the transaction priority number, with facilities pro- 
vided for the user to assign and change the priority of 
eecurrent. vuransact. Ol). 

During the setting up of a model, to be described 
later, the user is required to supply the number of trans- 
action parameters required. Each transaction generated 
in the model will then have the specified number of para- 
meters associated with it. These parameters are generally 
used to indicate uniqueness among the various transactions 
which may be in the system at any given time. For example, 
the transactions may represent ships in a berthing and un- 
loading system, and the parameters could specify variables 
such as ship size, type of cargo, and time to unload cargo. 

The transactions are maintained in two lists by the 
scheduling program. The current events list contains 
transactions whose next event times are less. than or 
equal to simulated time. The future events list contains 
transactions with next event times greater than simulated 
time. 

The current events list contains, primarily, trans- 
actions which are in a blocked state, i.e. they cannot 
execute the next activity due to the unavailability of 


some required entity. They will, however, be executed as 
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soon as possible. Also included in this list are trans- 
actions which will be moved at the given simulated time, 
including the current transaction. 

Transactions may enter the simulation in one of two 
ways; either generated automatically by the scheduling 
program or supplied by the user. To generate transactions 
automatically, the user must supply the value of two OSS 
variables, MEAN and MOD, with MOD < MEAN. The transactions 
will then be generated every MEAN + MOD time units, with 
equal probability given to each number in the range. Note 
that if MOD equals zero, transactions will be generated 
at a constant rate, every MEAN time units. All trans- 
actions generated automatically are assigned a priority 
of one, the highest allowed in the system. 

The second method permits the user to place any 
number of transactions on the future events list before 
a simulation run. However, the user is required to supply 
all the system and user perenecens for each transaction. 
This second method is useful in testing and debugging a 
model. The program also allows the two methods to be 
used simultaneously. 

After a transaction has been created it may remain 
in the system for an extended period of time during which 
its state may be active, blocked, suspended or terminated. 


Only one transaction, the current transaction, may be 
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active at any one time, and the parameters associated with 
a given transaction can only be referenced at this time. 

A later section will describe the relationship of the 
transaction states to the movement of transactions. Trans- 
actions which move into the final state, terminated, are 
not required in the system any longer, and are deleted 


from the system. 


4.2.3 Permanent Entities OSS provides the user 
with two permanent entities representing two pieces of 
physical equipment of the real world, facilities and 
Storages. A facility is a piece of equipment which can 
be engaged by only one transaction at a time. A storage 
is a piece of equipment which can be engaged by any number 
of transactions at one time up to some predefined limit. 
Although a facility can be occupied by only one unit per 
transaction, a transaction may enter any number of units 
into a storage, subject’ to.the above’ condition: 

If at any time transactions can not continue through 
the model, for: example, if a facility is occupied by 
another transaction, they will be placed in a first-in 
first-out queue advancing one at a time as the entity 
becomes available. Queue statistics may be obtained in 
such an instance by using the system supplied functions 
QUEUVEIN and QUEUEOUT. These provide statistics such as 


maximum transactions in the queue, total transactions in 
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the queue, number of transactions not delayed in the 
queue, and average time in queue. 

In addition statistics are automatically gathered 
for each facility and storage defined in the model when 
the functions FACILITYIN, FACILITYOUT, STORAGEIN, and 
STORAGEOUT are used. Facility statistics include total 
number of transactions that used the facility, average 
time transaction. in facility, and fraction of use of 
facility. Storage statistics include average, maximum 
and current contents of the storage. 

Facilities, storages and queues are identified in 
OSS by assigning to each entity type successive integers 
starting at one. 

In addition to the above entities the system allows 
the user to define any additional variables or data sets 


which may be required in a model. 


4.2.4 OSS Scheduling Program One consequence of 


using a material or flow based world view is the necessity 
of supplying, as part of the simulation system, a rather 
elaborate algorithm to schedule the occurrence of events. 
The complexity varies with how much of the scheduling is 
left to the user. For example, SIMSCRIPT contains a fairly 
Simple algorithm as the user is required to program the 
occurrence of each event by means of event notices. How- 


ever, languages such as GPSS, SOL and SIMULA remove the 


J 2 Saws 


Bd 


ent at beysleb- joi enotaadanant! to tedmun _ Su: 
.eveup at omits snorevs bra :euewe 


betedtsn Ulsottenoius eis aottetiste soivibbs at a 


. fade Lebom ont at peated sastote bas yttitost dose 02 
brs JMIBOAROTS?, , TUOLTIAIDAY , WIYTIALOAD snoktonst ent 
fatot ebuLont gottettste yiiitost .boay ots MWOSDAROTE 
gpereve .vttitos? edt fait ehnolttosensit to sedan 
to ga to qoigos1t brs ~ydtitost at nottosensit emit 
. frusen be em , S88T9V8 sbufont aotsetiate sgsi10de Sdemcareys 
= | ot : " s BBBIONS edt to atagtrios’ jnsTwWs bas 
mt pefttiaebl ors ‘eoustio bre eonstota ,estziitost 
sregesnt evteassoue saut ystine doss ot aninglees yd 220 
| ; 90 36 SabIsss 
4 ewolls moseve: edt eotitine evods sit o¢ notdthbs al 
atee sisb 10 2eidsiusv fenoterobe yas entteb ot 1S8u eds 


cee | 


. Lebom s ot bevtuper sd - dio tetw 


to eqnaupeenoo ‘end. agregar patente 220 ~#, + e cians 


__ythezesq, ont at way birtow. besesk woft pt) fetretsm 8 gntey 


: tetist Bs) genanye Kobsotinte maa to g1sq 48 ‘SS ee ae 


-edneve to sonsi1u990 odt efubsrios ot nistropls stexodete 
Ja aattubsrtoe ods ‘to towm so sew edinay wstxetgmos sat 


7 la yesakvon neve. to bon we 1 anew‘ does 16 
pitt: von ieee ne a as foue’ as 


2 
¢? 


47 


burden of the user scheduling all events by use of more 


elaborate programs. OSS uses the latter approach. 


4.2.4.1 Organization An important aspect of 

the OSS scheduling algorithm is a set of user supplied 
configuration parameters. These indicate the activity 
which is to be executed immediately after a given 
activity, and are supplied as the model is being assembled. 
In general only one configuration parameter is required 
for an activity, however, OSS will permit the user to 
specify a aac hienen and secondary parameter providing the 
activity contains a TRANSFER statement. This statement 
allows the user to change the path transactions may take 
through the system and will be described later. 

Operating in conjunction with the scheduling program 
are a number of control statements. They are included 
in the activities by the user to indicate to the schedul- 
ing program the status of the current transaction after 
an activity has been completed. They include the follow- 


ing statements. 


1. ADVANCE A The current transaction is to become 
suspended and placed on the future events list 
for A time units, after which the next activity 
is to be attempted. 


2. BLOCK The current transaction could not execute 
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the “nextPactivity!* Dt *is’ kept *fnethe ‘current 
events list in a blocked state until the blocking 
condition is removed. 

BL THOBD FA The current transaction could not execute 
tHe, NEXt Sactivity becausesof /aeb locking pcondition. 
lt "isoplaced on the future *events Olist for’ A 
time unitsVafter which *the ®activity is again 
attempted, 

4, ENDTRANS The current transaction is deleted 
fiwom Scheiecyrnent vevernts alistypsthe @count.io£h& com= 
pleteditransactlonsuls! increnented *by vone feandsthe 
simulation elimi tis. wecheckederonean vendtef tsimulation 
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Hoe to Operation. As previously indicated .in 
Section 4.2.2 each transaction has, as one of the.system 
Val bebles s oeb bene act Lon Stats sindicavom.~ vor Lrans— 
aCcuLoOvemin the Current events list the status indicator 
bas the svalue one or four iIndicatine the transaction was 


blocked the gbast time it was active or the transaction is 
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available to attempt the next activity, respectively. 
Note that if all status indicators in the current events 
list are equal to unity no further state changes can 
occur at this-simulated time. 

Assuming the simulated clock has just been updated 
the scheduling program will proceed as follows. Each 
transaction in the current events list, beginning with the 
first one in the highest priority class, will be selected 
in turn as the current transaction. Each transaction is 
moved through as many activities as possible until a point 
is reached at which it can cause no further state changes 
at the current simulated time. After no further trans- 
actions can become current at the simulated time, the 
scheduling program moves the transaction with the lowest 
EVERUBCINCm lc OMe vice rULULes LON UnemcCUrrenv eventes liste 
where it is inserted as the last. transaction within the 
correct priority class. The clock is updated to the time 
of this transaction and the process is repeated beginning 
with the first. transaction. with the highest, priority. 

In order to provide a more realistic dynamic model 
the above sequence is altered under certain conditions. 
For example, a transaction near the end of the current 
events list may release a facility which is blocking a 
transaction near the beginning of the current events list. 


The moving of transactions would be started over as soon 
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as the former transaction became suspended. Thus, events 
will occur at the earliest possible time and in the cor- 


rect sequence. 


Lee OSS ECONETOL Program The use of OSS is under 
control of a.main program through which the user specifies 
the section of OSS he wants to use. The control program 
will type ENTER COMMAND WORD and the user will respond 
by typing one of NEWMODEL, STORE, LOAD, RUN, NEWRUN, END. 
The control program then accepts the command word, vali- 
dates it, and if it is acceptable will branch to the 
appropriate routine. After completion of the specified 
routine control passes back to the control program which 
requests another command word. The command words and 


their function are described below. 


1. NEWMODEL As the command word indicates a new 
simulation model is to be set up. The program 
requests the names of all activities that are 
to be used in the model followed by the con- 
figuration parameters. The user is then required 
to supply the following: the number of trans- 
action priority classes, number of transaction 
parameters, number of queues, facilities and 
storages, and the capacity of each storage. If 


the user wants to place transactions on the 
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future events list they are entered at this 
point, otherwise he can specify the values of 
MEAN and MOD to generate transactions auto- 
matically. Note that both methods may be used 
in the same model. 

The program then requests the initial 
simulation time and the criteria for stopping 
the simulation run. The user can terminate the 
run in one of two ways: by specifying the length 
of time the model is to run or the number of 
transactions which are to be created and destroyed 
during the run. In the present program the user 
must supply values for both criteria; however, 
the effect of one can be eliminated by specifying 
a large value for it. The model is now assembled 
and control returns to the OSS control program 
which requests another command word. 

One other important point should be mentioned 
in connection with this command word. Since OSS 
‘does not operate within the APL system directly, 
the task of setting up a model is at the present 
time "semi-automatic", i.e. the user must add 
activity names to parts of the OSS. control pro- 
grams while under direct control of the APL 


system. The areas of interest are shown in 
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Figure 5.. The function BRANCH is used during 

a simulation run to call activities to perform 
the state changes. The array FCNS is used only 
during the assembly phase, i.e. when the NEWMODEL 
command word is specified. Note that.the array 
contains the same activity names that appear in 
BRANCH and in the same order, and that the names 
must. contain six characters including blanks if 
necessary. Then before a model can be assembled 
and run, the user must insert the names of all 
activities he will use in the function BRANCH 

and generate an array, FCWS, containing the same 
names that appear in BRANCH. 

STORE This command allows a user to store all 
the required information of a model, currently 
active in the system, so that it may be retrieved 
and run at some later time. The user is required 
to supply a set of unique names by which various 
components of the model are stored. 

LOAD This command allows. a previously stored 
model to be retrieved and made current. The user 
must supply a set of names, as called for, by 
which the model has been previously stored. 

RUN After a model has been established in memory 


(by either WEWMODEL or LOAD command words) this 
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command word will cause the model to be executed. 
It will continue until one of the limit criteria 
are exceeded, until the user presses the APL 
interrupt key, or until the APL interpreter is 
unable to continue because of an error. In the 
latter two cases control is taken over by the 

APL system and the user can examine any variables 
or arrays in the system, or make changes to any 
activities. The program can then be resumed at 
Lne pein of interruption. 

If the simulation runs to completion the 
control program prints the time at which the 
simulation ended and the number of transactions 
terminated. The user then has the option of 
obtaining a print-out of the statistics which 
were gathered during the run, after which another 
command word is requested. 

NEWRUN Before executing a model this command 
allows the user to reinitialize some of the 
model parameters. For example, transactions may 
be placed on the future events list, MEAW and 
MOD may be specified and the length and initial 
value of simulated time may be specified. 

END This command takes the user out of the OSS 


system and returns control. to the APL system. 
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4.2.6 OSS’ Functions “In addition to the transaction 
status statements described in Section 4.2.4.1 and the 
complete set of APL facilities, a number of. additional 
statements have been implemented as APL functions to 
perform many operations required in simulation models. 


A description” ofeach, and their syntax, follows. 


1. .+(FACILITYIN A)/LABEL This is used when an 
attempt is made to seize facility A. If the 
facility 1s occupied the current tee enion 
is placed in a blocked state and the next 
statement in the activity to be executed is 
LABEL, where LABEL is a valid APL statement 
label. Otherwise the facility becomes, occupied 
and the next sequential statement in the activity 
will be executed next. An alternate form, 
R+FACILITYIN A, where R is a dummy variable, 
may also be used. However, if the facility is 
occupied the next sequential statement will be 
executed next. 

2. FACILITYOUT A Facility A will become available 
to another transaction. Corresponding statistics 
will be updated. 

3. +(A STORAGEIN B)/LABEL This statement will add 
B units to the current contents of storage A 


unless the storage capacity will be exceeded. 
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If exceeded the current transaction is placed 

in a blocked state and the next executed state- 
ment is LABEL. Otherwise the storage is entered 
and the next sequential statement is executed 
next ., The alternate form, R+A STORACEIN B 

may also be used with the same result as described 
imekbmeabove:. 

A STORAGEOUT B B units are removed from storage 
A and the storage statistics are updated. 

A QUEUEIN BB units are added to the current 

and total contents of queue A. Queue statistics 
are updated. The transactions are treated on a 
first-in first-out basis. 

A QUEVEOUT B B units are removed from queue A 
and the relevant queue statistics are updated. 

A -ASGN -B ‘The value of parameter :A..of the current 
transaction becomes B, where A is less than or 
equal to the number of parameters specified in 
sevting sup tthe emodel, tandsB is,,any «integer 

within the capacity of the machine. 

PARAMETER As This; twill retrieve athe value of 

the ee parameter, not exceeding the specified 
number, of the current transaction. 

PRIORITY A. This will assign the current trans— 


action ia priority A. 
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TRANSFER EXPRESSION. As described in Section 
4.2.4.1 this statement is used in conjunction 
with a primary and secondary configuration 
parameter, and allows the programmer to change 
the path taken by various transactions through 
the model. The value of EXPRESSION determines 
whether the primary or secondary configuration 
parameter will be used as the next sequential 
activity seone [for stheseprimarysand azero wfor hthe 
secondary. configuration. Thus, EXPRESSION can be 
any valid slogical APL expression. For example, 
TRANSPER (QCONT «1)"s65, will cause the activity 
given by. the primary configuration parameter to 
be executed next if the contents of queue 1 is 
less than six, otherwise the next activity will 
be given by the secondary parameter. 

FACIDLLY. AncThesresult elsisthexscurrent tstate: of. 
haciiity Ase letter eccupled,, 0 otherwise. 

STCONT cA The result “is: the currenti«contentsyvof 
storage A. 

CONT MAM =InemresiultelSa vie current! Conventsr ot, 
queue A. 

A RAND B This statement generates a number in 
phe range A.+t 8,(B s A), with equal probability 


given to each number in the: range. 
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4,3 An Example 
To illustrate the use of OSS a model will be con- 


structed illustrating the operation of. a simple telephone 
exchange. [elhevisystemrcons¥s ts scofea munber eof telephones, 
EGaGm@.connecved to the exchange by its own) line. © The 
exchange contains a number of cross-links which can provide 
a connection between any two telephone lines entering the 
exchange, with only one connection being allowed to any 
SLVenmelineeatteastime.e subs s assumed that calls which can=— 
not be made at a designated time will be attempted until 
finally made. For simplicity assume the system contains 


EO telephones “and 5 cross—links. 


4.3.1 Developing a Model To develop a model of the 


system the following are necessary: 


1. Entities: to) indicate. ‘the *begin' and end time of 
each call and the two telephones involved in the 
Ca li 

2. Entities to indicate the state of the telephones. 


3. Entities to indicate the state of the cross=<links. 


The first entities can be recognized as temporary and 
will be represented as a transaction in the system, with 
the telephone indicators as transaction parameters. 

Since the telephones and cross-links can be in only 


one of two states they can be represented as facilities. 
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ple) 


Then before a call can be placed a check must be made if 
the telephones are available and if a cross-link is 
available. Queue statistics will also” be “obtained on 

those calls which cannot be placed at the earliest possible 
time. 


The entities to be used are defined as follows: 


Parameter 1 Telephone number sending call 
Parameter 2 Telephone number receiving call 
Parameter 3 Length of telephone call 


Faciiicy ~l- Telephone number 2 


Facility 10 Telephone number 10 


Bacidity—ll Cross—Link number 1 


Facility PE5emtross=Link umber?! 5:. 


Telephone calls are placed every 10 + 5 minutes, 
their length being 3 + 2 minutes. One unit of simulated 


time will represent one minute of real time. 


4.3.2 OSS Programs Excluding the activity to gener- 
ate the transactions; three activities are required as 
shown in Appendix II. The following is a description of 


the function of each activity. 
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2. Assign the telephone numbers used 


in this call to parameters one and 


two by use of the random number generator. 


Check to ensure different telephones are 
involved in the call, otherwise go to 
line 2 and reassign the telephone number 
receiving ltheicaLlt 

Assign parameter three the length of 

the call chosen from the range 6 + 2. 
Phacesa@this call into a queue “of calls 
to be made causing the queue statistics 


to’ be updated. 


Checkssthe availability ofboth. tele- 
phones required to initiate this call. 
Lim Ones Ls netmavailab le @lhe cal las 
placed. in a blocked state at line 10. 
This line generates a vector ZINK con- 
taining the facility numbers correspond- 
ing to available cross—links.. if-no 
Gross—linksmaresavar lable tnescalli@is 
placed in a blocked state. 

Removes the call from the waiting queue. 
5. wets the facilities corresponding to 


the required telephones to occupied. 
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Line 6. The first available corss-link is used 
to place the call. 

Line 7. The number of the cross-link used is 
placed in transaction parameter four. 

Line 8. The transaction is suspended an amount 
of time equal to the length of the call. 

tine 9. “Exit from the’ function, 

Pinew Ove lacess cad I seinealp locked oteale. 

3.  ENDCAL 

Line 1 & 2. Releases the telephones used by 
CAs Se calla 

Line 3. Releases the cross-link used by the 
Carlie 

Line 4. Deletes the current transaction from 


the system. 


4.3.3 Results After the activities described in 
the previous section were written and entered into the 
APL system, OSS was used to assemble them and run. the 
resulting model. Appendix II contains a complete listing 
of the setting-up, running, and statistical output obtained 


from the model. 
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CHAPTER V 


SOME APPLICATIONS OF OSS 


In this chapter we will consider three systems 
which have been simulated using OSS. A description of 
each system will be given followed by a discussion of 
the simulation model developed for each system. Through- 
out this chapter STU will be used as an abbreviation for 
Simulation time units and the real time corresponding to 


one STU will be specified for each problem. 


5.1 Simulation of a Three Machine Assembly System 


This problem is the same as that presented in 
Section 3.4 where it was simulated using the FORTRAN 
program. It has also been simulated using OSS to illus- 
trate the differences in the two approaches from a pro- 
gramming point of view. Refer to Section 3.4.1 for a 
description of the problem and to. Appendix III for a 
listing of the OSs activities used to simulated the 
proulem, In this example one STU will represent one 
MiInuce Or real time. 

Since the system contains components produced by 
machines A and B as temporary entities, these will be 
represented as transactions in the model. However, 


since they are produced and used in pairs, each trans- 
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action will represent two components, rather than one. 
Random numbers are used to determine if a component is 

to be rejected, i.e. an integer number in the range [1,10] 
is generated and if equal to 10 a component is rejected. 
Then, because one or both of the components of a trans- 
action may be rejected, we must ensure that some of the . 
transactions are terminated when a component is deleted. 


Three cases may be recognized. 


(a) If both components are rejected the transaction 
corresponding to them is terminated. 

(b) If neither component is rejected, the trans- 
aculOnei senoOuscermineved. 

(c) If one component is rejected, terminate the 
transaction only if there is an even number of 
components in the storage bin, denoted in the 


model by STORAGE 1. 


The number of components to be entered into the 
store by a transaction is indicated in transaction para- 
meter 1. .As shown in the first’ activity, REJECT,«the 
first statement assigns parameter 1 a value of two. 
Subsequent statements,in the, activity decrease the 
parameters bypone| fornevery {component rejected.» The 
second half of the activity contains statements to deter- 


mine the number of components to be entered in the storage 
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bin, to add that number, and determine whether the trans- 
action should be terminated or not, using the above 
mentioned tests. 

After a transaction has completed REJECT and has 
not been terminated it immediately attempts to seize 
FACILITY 1, representing machine C... This action is 
represented as the first. statement. in activity MACHNC. 
LS the, facility. is,.occupied,.the .transaction 1s, placed 
in a blocked state and the activity is immediately 
exited. If the-facility is net eccupied it is' seized 
by the current transaction, two components are removed 
from the storage bin, and the transaction is suspended 
for 15 STU, representing the processing time of machine C. 

After 15 STU the suspended transaction is placed 
back on the current events list from the future events 
List,-and. when <it..becomes -current,..the last activity 
ENDTRN is executed. The state changes caused by this 
activity include incrementing the count of completed 
products, represented as STORAGE 2, by one and releasing 
facility. -one.¢ The jtransaction is then deleted from the 
model. 

The activities were assembled as shown in Appendix 
III. Note that in this and other problems which use the 
automatic generation of transactions the first activity 


specified is GENRAT, although no activity of that name 
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is actually executed. It is, however, necessary to allow 
the control program to set up the correct activity sequence, 
as well as providing the user with a listing indicating 
transactions were generated automatically. 

The statistics obtained indicate 19 products were 
completed during a run simulating an eight-hour shift. 
During this time the maximum number of components in the 
storage bin was three and machine C was in operation 


SOO fe of the time, 


5.2 Simulation .of a.Real-Time Data Processing System 


The system being simulated in this section consists 
of a computer with three disk drives sharing a channel 
to the CPU, and a high speed on-line terminal. The computer 
receives messages from the terminal, and a search is made 
en one of the disk files to retrieve a record which is sub- 
sequently used by the computer for processing. Figure 6 
illustrates the sequence of actions that are required in 
the model. In this example one STU is equivalent to one 
millisecond of meal, time. 

Messages are sent by the terminal to the computer 
every 50 + 25°STU at which time the’ CPU-requires 1 STU 
to place the message in memory. Each message enters a 
queue since the disk file can search for only one record 


at a time. When the appropriate disk becomes available 
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the search request corresponding to the message at the 
top of the queue takes control of that disk. Before the 
correct record can be retrieved the following two cone 


ditions must be met. 


(a) The disk arm must be positioned over the correct 
CrackteOustier disk. se liicmheduiresm! cOmrsoOr STU. 
(b) The correct record must be under the read head. 


Tiiserequiresmebr 2 25oTU. 


After satisfying these conditions the record may be 
read into the CPU if the channel is not being used by 
one of the other disks. If not available the attempt 
to read the record must be suspended until one complete 
revoluclonsots une diskew or 50, pTUy = ihis dserepeated until 
the channel is available. When the channel becomes avail- 
able ath Sse seized@for=> 9sTU tov¥pertorm the read. The 
records cane then ber processed whensthe CPU is available, 
after which the message is removed from the computer. 

The temporary entity of OSS, the transaction, is 
used to represent the messages which are sent to the 
computer from the terminal. Each transaction has 
associated with it two parameters: the first indicating 
the length of the message, and the second indicates which 
of the three disks contains the record of interest. The 


transaction then, moving through the model, represents 
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the action of the equipment through time in response to 
a message being entered into the system. The system can 
be used to study the effect of different message intere 
arrival times on the time required to process a given 
message due to the single channel. 

The system equipment has been represented as the 


OSS entities listed as follows: 


Racw Uy @lareDa Sle lUrister 
2 MKC W Ee Sule aly, 2 
See CUskeuniees 
4 Channel 
a CP 
Storage 1 CPU storage for messages, capacity 


1500 storage units. 


In the example shown transactions are created every 
50 + 25 STU. The message length parameter is chosen from 
Lhe rectangular, distribuciomsl5iiel0. Associated with 
each disk is a queue, e.g. queve 1 with disk 1, which will 
be used to obtain statistics which can be used to judge 
the system performance. Following is a description of 


the activities; listedsin Appendix LIL; used insethe model. 


1. MSGIN Each new transaction can enter the system 
only when the CPU is available. When available 


this activity assigns the two transaction para- 
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meters and suspends the transaction for 1 STU. 
FACOUT This activity releases the CPU facility 
to another transaction. 

OLOREN ialfiealhethe CPU sittorage ‘is’ occupied sby 
previous messages the transaction will be blocked 
until some storage is released. The activity 
puts the transaction in the queue indicated in 
parameter two, after being entered into the CPU 
storage. 

Qifst welt thevdisk Taciiivysis,occupledsthe 
transaction remains in the queue. Otherwise 

it is removed from the queue and suspended for 
120 + 80 STU, representing the time required to 
Post tilon® (ne) read neadeoversthe correctstrack. 
ADVANC The transaction is again suspended, but 
Forse) taco) olU. Corresponding  cogwaltinge= for 
the disk to rotate and position the correct 
record under the read head. 

READIG® After positioning the-recerdrin the 
previous two activities, this activity first 
checks for the availability of the channel. 

If not available the transaction is suspended 
for 50 STU when the activity is again attempted 
as indicated by the line HOLD 50. Otherwise the 
channel is occupied for 5 STU representing the 
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( SENDTSHeeThis activity releases: the disk and 
Channel facilities to other transactions. 

8, PROMES A check is then made on the CPU facility; 
if available the CPU is seized for a calculated 
time as the record is processed, otherwise the 
transaction is blocked. 

9. ENDMES The CPU facility is released, the message 
is removed from CPU storage, and the transaction 


is terminated. 


The results obtained for this example are shown in 
Appendix III. The run was completed after 2995 STU during 


which 47 messages had been completed. 


5.3 MoLrmulatilon of “a small Branch Library 


This example is taken from Caras (1968) and is out- 
Paces eire Vee lo as divided mainve two basic sections, 
books and documents. 

People enter the library every 5 + 2 STU and leave 
if there are more than eight persons waiting at the infor- 
mation desk. The people at the information desk are 
served on a first-come first-serve basis, each requiring 
8 + 3 STU, to be directed to either the books or documents 
section. 30% of the customers are sent to the books 


librarian, the remainder to the documents librarian. 
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There is only one books librarian, and if there are 
two people waiting for service the customers diverted to 
this section will leave. Otherwise they form a queue 
requiring 25 + 5 STU for service. After being served 40% 
leave the system immediately while the remainder attempt 
to use a duplicating machine. If it is in use the customers 
leave the system, otherwise they use the machine for 45 + 10 
STU, after which they leave the system. 

Customers diverted to the documents section divide 
evenly between two documents librarians forming queues at 
each. The first librarian requires 18 + 3 STU to serve 
customers, the second 22 + 4 STU. After leaving the 
documents librarian three-fourths of the customers leave 
the system. The remainder attempt to use a film reader 
which is in a small room with a capacity of three persons. 
If the room is full the customer leaves, otherwise he joins 
a queue for use of the reader. The person using the reader 
requires 35 + 5 STU after which he leaves and the first 
person in the queue takes over. 

Since we are obviously interested in the movement 
Of people through the system each transaction will repre— 
sent one person in the system. The various entities 


required Inevuhe moder are derimedwas follows” 


SUSU 8 in6% ‘yout eehoy esa ee 
You favrse goisd -totTA -sptvass aot Ute @ 2 8S 
sonsiss tebritsmst eid obeige ‘yletstbemnt neteye ‘odd 

etemod ao edt. eau nt af ¢b tL entitssm aitttsotiquh s Bae, - 
OL + tee ‘to? eftdoem eft sau youd seiwrerito ,mSsteye ent svsel a 
| smeteye “odd avast yedd dotdw setts OTe 7 

spivis nottose edmemuoob ent od pettev tb atomoten) | bea 
% seve jatmso ensirerdtl atnemuso0b ows: “geewied vines 7 
avise: ot- ute € + gr gor huport netistd tt sank? ont? 988 i 
~sid gatvaet tests , .UT2 ts ss brossa oft < STOMOTEHD ; 
ovbet ergata eit to. eridcuot-oontd asiretdtt esreew ood. | 


xebeer. mitt a seu 02, tqmedts sobatatey oft nodeye at 


. 
enoRtsg otis ‘YO “Yttosass 8 ailw moot {fsma & ato dtd 
entot ‘ed: petwrenito ,noves! aemod ans ott ‘fiw at moor askin i 
tebse1 odt gntey moeteq ont, 19b89% adit 38 . 
2 q ath x02 a ay! 
tent? and bas eaves! ed doteiy setts pre-e + ee a 

stevie eoxst eee 
sitentsvor ont nt  batesroent vlenotvde 9) S18 wt: 
i a £Ltw aoptossnesd dose mopees re 

) ha ya 


. 


fe 


Facility 1 Information. desk 
2 Books librarian 
Duplicating machine 


First documents librarian 


U1 


Second documents librarian 
6 Film reader 


Storage 1 Film reader room, capacity 3 persons. 


In addition each of the above facilities, except 
facility 3, may have lines of people queued up and waiting 
for service. To obtain information on the waiting times 
and number of persons waiting for each facility, a queue 
Wille besassoclated with each Tacitityv,, =the statistics 
obtained can be used to analyze the efficiency of each 
section and to pinpoint any bottlenecks in the system. 

Following is a econ nnten of the activities required 


iethis model, sas =listed in Appendix Ii1. 


1. ENTER. The-contents:of-the information desk 
queue is checked; if more than 8 people are 
waiting for service the customer: leaves and the 
transaction is terminated. Otherwise the trans- 
achmuonea.s placedsain) queues ly. 

2, WNFORMe This activity represents the action of 
theyinformatvions desk, eThej first dine.checks the 


status.of.the facility. If available the trans- 
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action is removed from the queue and is suspended 
8 + 3 STU representing the time required by this 
customer at the information desk. Otherwise the 
transaction remains in the queue. 

SPLIT This activity provides the necessary 
statements to accomplish three tasks. First a 
random number is used in conjunction with a 
TRANSFER statement to send 30% of the trans- 
actions to the books section and the remainder 
toathecdoecuments (sections asSecond, the trans— 
actions routed to the books section check the 
contents of the waiting line at the books 
librarian. pifinvheresare nobmore rthan twospeople 
waiting for service the transaction is placed in 
Therqueue. Sotvnerwisegthesuransactioneirsaterminated 
indicating thercustomers leaving the system. Third, 
transactions routed to the documents section have 
a parameter assigned indicating which of the docu- 
ment librartanssis to service the transactions. 
BOOcKL I Gylinthe-bpook*itbrarian Iskavatiable the 
transaction is’ removed: from the corresponding 
queue and suspended for 25 + 5 STU, the service 
Bimetor this Librarian. 

DUPMAC This activity terminates 40% of all 


transactions entering it, representing those 
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people who had no use of the duplicating machine. 
The remainder will also leave if the duplicating 
machine is not available. If the machine is 
available it will be engaged for 45 + 10 STU. 
ENDTR2 Transactions leaving the duplicating 
machine release the facility and are terminated 
Serie Wits: Tarot 0 ya: 

DOCLIB This is the path the majority of trans- 
actions will take after activity SPLIT. Each 
transaction carries in parameter 1 the queue in 
Whichyit was placed ins SPLir, with queues: 3 and 
4 corresponding to the first and second documents 
iibprarian, presvec ul yely she tirsy tine or this 
aeviv teyecheeks sUhewavawlap i! Loy sOneaLnewredulred 
ibrariant and it avallaple the transaction is 
Suspended a period of time selected from the 
distributaon corresponding to the required 
allshachouliche 

FILMQU After releasing the document librarian 
facility three-fourths of the transactions are 
terminaved. The remainder check the number of 
people waiting to use the film reader, and if 
greater than two, they also are terminated. 
Otherwise the transaction enters the film reader 


queue and storage entities. 
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9. READER A transaction finding the film reader 
(facility 6) available will be removed from the 
queue and suspended for 35 + 5 STU representing 
the use of the reader. 

10. ENWDTR1 The transaction leaving the film reader 
is removed from the storage and reader facility, 


and is terminated. 


Figure 8 shows the manner in which the activities 
are joined together to form the model. The integers 
associated with each activity are those used as the 
configuration parameters when the model was being assembled. 
Note? thar *ror™“activity SPLIT the primary and "secondary con- 
figuration parameters were specified as 5 and 8, respectively. 

After 100 transactions had been generated and termin- 
ated, the results shown in Appendix IIIT were obtained. It 
is apparent that in this: hypothetical system facility 1, 
e€orrespondine vo the information desk.,jcreatves a bottleneck, 
since on the average people had to wait in line 67.81 time 
units before being served. Thus, better service at the 
information desk is necessary. On the other hand only 7 
people used the film reader indicating another is not 


required. 
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CHAPTER VI 


CONCLUSIONS AND SUGGESTIONS FOR FUTURE RESEARCH 


OR eerconclusions 

The following conclusions are based on a comparison 
of the two simulation programs presented in this thesis 
with the simulation language requirements specified in 
sections 2.1, 2.2, and 4.1.2, and on experience: gained 


in using the programs. 


6. . eeeeORUR ANRC umu le tL Bon Program Although the 


Oroeram providedea majomkty Ol the eeecuimements of. ra 
Simulation language the following observations may be 


made. 


1. The arrangement of all entities into one large 
array, With the same name for every type of 
envoy. LSeawkward TOUUSeetrom a users’ point 
of view. A better approach would be to provide 
aenumbereotmoypes Of ENUlUvLeS, se .e. Lacilipies 
Bho vOcCases wall eCoCiereLlerenced. by sae name 
Gescriptave oL thav vype. 

2. Although a machine-based simulation language 
Pequires sno alsorlunmevo, schedule the activities 
some mechanism must be provided to advance simu- 


lated time. The fixed time increment used in this 


e 


HORARCIA AAUTUE AOT 


Seatac’ £ a0 beead Ts Co antwol io" adit 
etesdd eidt mi patebeaiie -ehstaotg notte tome owt adit 
nt betttsoge _ainsmenLupsty sBei/anst nottalumbe ead ote 4 7 
heniten sors tiagKe” ‘fo ‘faei eS. Tit Sra, S8s5) Ye aii 
-poneeigo%a al anteu ou " 


ad 


y x ' ; 


vy 


= ae 


ehh eno odat pelsiene {Is to shombansta sot - 

to say? ytevs ot smart emse.. od dtiw vera 
| intog, Sane 6 mont Sau ot puswiws ar. uotine: are 
es oe es od ‘bisow dosotrqas netted o% orate 3 | a 
| abteEEos? aes (Setdtine st noay? to. bcacaersali : 
om 8 Nd, Peousiotos oss aoew cabaarote bos 
7 ae, adie denis” 0 ee 5 
¥ WE ee posites posads oe ‘Bt 
i _ sebdivison ent ‘tiboinn 20 9 ate 


Se 


% 


ihe 


program is usually wasteful of computer time ag 
activities will be attempted at times when no 
Steve changes twiltitoccevume liThifs tis: airesul tof 
making the time increment small enough to ensure 
that ale events will occur at their proper time. 
The use of special time entities such as described 
inwPeculoOns2.5.1l.1 would Overcome this disadvan— 
tage. 

One essential feature of a simulation language, 
which this one is lacking, is the automatic genera- 
C2LOn sor "stvatrstl cal wore ouner performance: data as 
UReCeMOdel Ss Dene run selhted See COnSequence OL 
MS) Weiqhwaien’ Gye eiee Phas. ivsiahe nel@gisel ake Wey lefeniicrs “sikaKels 
Chen user is rece ro deiiner an snvlty TO represent 
any required variable. It may be noted that a 
relatively large amount of data is obtained from 
the OSS simulation compared to that obtained from 
the FORTRAN simulation for approximately the same 
amount of programming effort. 

Although a simulation language may provide all 

the features’ as indicated in Section 2.1 ‘and 

One peulie TUtIkGy wor the: Kanguege wWependsy on diow 
Gasunky Lilet VETIMTeavures Cai be used VOM cOonstruct 
Geeaimulacion model andton tow tréediistic«the etinal 
model may be. Thus, although the machine-based 


world view of this program permits the construction 
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of realistic models, many of the necessary 
features of a simulation language can only be 
obtained after considerable effort on the part 


of the user. 


Sei Mies 


1. The OSS program provided all the simulation 
language features, including those of on-line 
systems, with the exception of a compiled mode 
of execution as this feature is snot currently 
implemented in the APL system. 

2. OSS proved easier to use than the FORTRAN program 
because: 

(a) it allowed the separation and referencing 
of entities by meaningful class names, 

(b) numerous statements were provided to perform 
operations on the entities, and 

(ec) the simulation models were written, debugged, 


and assembled in a time-sharing environment. 


3. Aseseensinethe,applicationsof OSS, the use; of 
the event as the basic descriptive unit results 
in some activities that contain very few state- 
ments. This could be improved by incorporating 
thersprocess as the Géscriptive unit. However, 
this implies modifying the APL system to permit 


"narallel" execution of the activities, in the 
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Same manner modifications permitted SIMULA to 

be incorporated into ALGOL. 

The powerful vector and array operators make 

APL an excellent language in which to write the 
OSS programs and the simulation models. 

A detailed comparison of the two simulation 
programs presented in- this thesis would be very 
difficult since they are written in different 
languages and operate in different environments. 
However, the following observations may be made. 
Of the two world views, material-based of OSS and 
machine-based of the FORTRAN program, the former 
was the easier in which to develop and write a 
simulation model. This can partly be attributed 
to the fact that ane systems are more amenable 
to simulation by a material-based language than 
others. However, the deciding factor was the 
numerous utility functions, built-in statistics 
collecting functions and division of entities 
into classes available in OSS. In addition the 
ability to create and destroy temporary entities 


is a very useful feature in a simulation language. 
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6.2 Suggestions for Future Research 
1. With the growing availability of on-line systems 


and the increasing complexity of simulation models, 
the need for an on-line simulation language will 
increase. Perhaps an attempt could be made to 
incorporate the OSS features directly into the 

APL system or as a subset of the language. 

2. As mentioned in Section 6.1.2 the process repre- 
sents a more desirable descriptive unit than the 
event currently used in OSS. Any plan to incorpor- 
ate the OSS features into APL or any on-line lan- 
guage should include changing the descriptive unit 
to the process. 

3. A relatively new piece of hardware used in time- 
sharing systems is the graphic display or CRT. 
Investigation could determine the usefulness of 
this input-output facility operating in conjunc- 
tion with an on-line simulation language. The 
CRT could be applied to problems such as drawing 
flow diagrams, building activities, input of 
parameters, and output of simulation results either 


during or after simulation runs. 
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APPENDIX I 


LISTINGS OF FORTRAN SIMULATION PROGRAMS 


This section contains a listing of the FORTRAN simu- 
lation programs, consisting of the mainline program and 
the subroutine INIT. Following these are listings of the 
subroutines described and used in the example of Section 
3.4 and a listing of the output obtained when the example 


problem was run. 
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CONTROL PROGRAMS 


VENERAL SIMULATIUN PRUGKAM 


NuUSty dS DNPOT UNIT. DF /S GR 

MON; 1S UYLTPUT. UNIT, Uscv TU MuUNITUR UdteRS TiNPUT 
NOUdw lS. CUTP UT) POe. USER 

NUAT 1S CATA iNPUT UNIT DURING SIMULATIUN RUN 
MAG IS MALTAPE UNIT 

ACTIVITY NAMES AKE 4 ALPHAMEKIC CHARACTERS STOREU 
AS NAMOP Ti) sl=iylo 

Hepsi TH AGI Vdaiye UF ANY MUveclL 1S OF KINDIT) AND 
NUMBER NUMBIT) 

PEACH MUUEL 15 STOREG UN FOUR ALJZACENT RELOKUS 
FINAL. RECUR DUS IS. NG eg PERiNbie AATINVAT I cS, jNup 
VERINGY hUUrLS, AND NAMES UF DEFINED ACTIViITI£ES 
NUPS=NUMBeER UF LEFTINED ALTIVITIES 

NCLUN=NUMocnK UF DEFINED MUUDLLS 


VIMENSIUN STATEMENTS 


DIMENSTUN NAMUPL1 LO) »yRINDI 144) ,NUMb1144) pRPAKRI LOU), 
LNKPAR( 16) pNPPAK( LO) psTATE( 25540) pPPARIiV00) 
DIMENsIUN LlwCRD(2) 


COMMUN/ BLKA/ STAT, PPARs KPANg NAMUP gy NAPAR y WPPARs ALND yNUMD 
CUMMUN/ BLKD/NDLUE NSTATy NAP yNPP 2S KP sKP Ps INDEX 
CumMmMON/ BLKU/S NULN INUPS »>MAGs NUSce yMUN rNUJUT yNDAT y TIMe 


a) SPECIFY UNIT NUMBERS ANL CUMMANU CUVES 


NUSE=5 

MUN=o 

NUUT=o0 

NUAT=5 

MAG=3 

DATA KUDE]L, NODE, KOLES ,KCDES,KUDES-, KUDES pKUDE/y Ju *RUN ', 
PACA 3 oP SUR yp 4 NERA sas Nie Maclin * ys. * ae 7, 
BDALOPENDVAL/ C2181)! o/s NOS TEND 4/7 


StdewNUoLU = —~lL 10 \INUVICATE NU MODEL IN, STURAGE 


NUSTU=-1 
NEG=-1 


Sha NUUN=—-1l TO INDICATE MAGTAPE HEAD MUST BE SET 
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NCUN=-1 
WRITE(NUUT,101) 
101 FURMAT (4,145) 
C 
& (2) REQUEST A B5ASIC CUMMAND 
c 
200 wKITeE(NUUT,210) 
210 FURKMAT (* C) ENTER RUN/CALL/S STURC/NEWALT/NECWMUD/® 
L'NcwTAPE/eND/INIT/*) 
READANUSEs 220). 1WORD 
226 FURMAT(2A4) 
WRITEC(MUN, 230) IwukD 
230 FURMAT (# UJ 4,2A4) 
ITEM = iwURvi1)d 
DEL IEM ER eNDde STOP 
Prldsiiet bUSCDESG) GO TO 6uU00 


€ 

FF PUsiTUON, REAG HEAD AT END UF MAGTAPE ITF NUT DUNE 
C PREVICUS io¥ 
* 
2 


‘abe fa TAANGUN, Nioeilisid) GU) aoe bo 
KEWIND MAG 

240 REAL (MAG) 1 
[TfFiLNE.(-1))60 TU. 240 

24D BALK SPACE MAG 
READ{MAG) WEG,yNOPS»)NCUN, NAMUP 


5 

. GRANCH 7.Cr APPROPRIATE RUUTING DEILEKRMINeED BY BASIC 
ve CUMMAND 
C 
Z 


250 IFCIT EME ACD EDO GOydIG LOGO 
lFliievteenm KGOe 2d GO Tu 2000 
LEAD LEP Eee A UUES) GU 1G) S000 
lFodviiche ree nhUvcec4) GU TL 4000 
PAT Mec etewcody GU lie > OGG 
Liisi T Gite Gees ALA, CALL INIT 
GU al Us 20:0 


G 
atc ti) USER HAS, ENTERED NEWT APE 
u PLAT SINGLE -LUGICAL KEGURD UN MAGTAPE 
C (-1,0;0,16 GLANK WURDS) 
L 
eCOG REwIND HAG 
NUPS=0 
NCUN=0 


DO 6010 J=l;1lo 

601G NAMUP( 5) =KUUc7? 
WKITE{ MAG) NeGyNUPSyNUCUNyNAMUP 
Gis gett «20,0 
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(4) US cn tino CNICKRED NEWALT 
Reeve So haMise) Clits THEMs STURGES UN LAST 
MeatG iG icyie) We vate sy ere 


SG COG Gy 


UCGO BALKSPALE MAG 

WRITE (NLUT;,4010) 
SsUbO SFoRMAI(' COC) ENTER. S“CHARACTER NAMES UF NEW ACTIVilledS*) 

wWRITcCUNGUT,4020) 
4020 FURMAT(! Cin BACH FOLLOWED Bye SINGLE BLANK*) 

ITEM=NUPS+4#1 

READ(NUSE,4030) (NAMUP{J),J=ITEM, 16) 
4030 FURMAT(L0(A4,14)) 

WRITE(MUN, 4040) (NAMOP(J),3=iTEM,16) 
4040 FURMAT(4 U) *FyLl0O0CA4,1xK))} 

WRITeE(NGUT 4045) 
AO45> FURMATHZ, C) SEFUORE DEFINING THEIR CUNFIOUKATIUN, COMPILE!) 

DO 4070 J=ITEM,16 

LE CNAMNGP (J). ESe RUDE) GU 10 40 7¢2 

NUPS=NUPSt1 

WRITECNGUT,4U600) JyNAMUP( J) 
4000 FURMAT(! ty) SUCKOS Nee PN si ce ® FUR ACTIViTY © 3A4) 
A070) CONTINUE 
4072 wRITECMAG) NtGy,yNUP Sy, NCU 7NAMULP 
GU. fU_ 200 


C 

C (5) 2SUSER GAS, ENT ER EDS NEWMUD 

< RERUESTUNS AIT Sono. GF STATE ENTDITES 

c Rewies. NNDLUK = NUe UF ENTITY 6GlLocKks 

e REWUEST NiDT = Nu. wr ALTIVITIES 

e REGUS TSCRUENCE WReUPERATIUNGS 

C SE ie dime l=_ 050 

L 

5000 wRITE(NUCUT,9C610) 

S0n0: FORMAT <4 CY) Shia dae NUMBERS Crass Aes che TT LES aANISs 45 


iS eecCK om Net CRMAL et2i2a 4 
REAVINUSEs5020) NSTaT,NbLUK 
5020 “FURKAl fiz) 
WRITE( MON, 50302 NSTAT,NDLUK 
5050  RURMAT<? OO) Sf oxt2s 
TIiMERHRaaC.0 
wRITEINGUT 50553) 
5055. i.ewewaAd af SO Seen Oe oe ALILVITIES iN MUDEL EN, 
lene Caveats 115743 
REABCNUSE, 5026) 4NTOF 
526  EORRAT IIS? 
WRIPEL MON, 5037) NIG 
S@3% FORMATAS Be) *,i15) 
wRITe({NGUT,5040) 
5040 FURMAT(! OO) Ru Se SseRVUENGe Ur ACHAVITITES "2 


3100 


5164 
51lod 


5los 
Si79 


DLeU 
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READ(NUSE,5050) ((KIND( J) »NUMG(U)) pJ=l,NTUT) 
PORMATHA SLAG 11, 1X)3 
WRITECMGNs5055)- GUCRKIND( 2) NUMe i) ) aL, NTUT) 
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APPENDIX II 


LISTINGS OF OSS PROGRAMS 


This section contains listings of all the APL programs 
required by OSS. They have been divided into the following 


sections. 


1. Control Programs The function SIMDRIV accepts 
the command words and branches to the appropriate 
routine... The other functions in this section are 
called directly or indirectly by SIMDRIV_ to per- 
form the functions indicated by the command words. 

2. Scheduling Programs The programs of this section 
are used only when a simulation model is being run. 

See eansaction status Programs These programs set 
the global variable ZSI which indicates the state 
Oiernececurrenvetransaction. 

4, Utility Programs These programs are the OSS 


functions described in Section 4.2.6. 


Following the listings of these programs is a listing 


of the activities and results for the example of Section 4.3. 
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obtained when the systems described in Chapter V were 
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VREJECT(COIV 


REJECT 

1AASCHT 2 
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1 ASGN( PARAMETER 1)-1 
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1 ASGN( PARAMETER 1)-1 
>+(O0=PARAMETER 1)/12 
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COCs, CON) so eh CSTCONr 1)22)/11 
R<+1 STORAGEIN PARAMETER 1 
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R+1 STORAGEIN 1 

ENDTRANS 


VMACHNCLOIV 
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+(FACILITYIN 1)/0 
1 STORAGEOUT 2 
ADVANCE 15 


VENDTRNCLOIV 


ENDTRN 

R<«2 STORAGEIN 1 
FPACILITYOUT 1 
ENDTRANS 
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ENTER COMMAND WORD 
NEWMODEL 


LYPE NAMES OF ACTIVITIES TO BE USED IN THIS MODEL 
HACH NAME CONTAINS 6 CHARACTERS SEPARATED BY A BLANK 


GENRAT REJECT MACHNC ENDTRN 


TYPE ONE CONFIGURATION PARAMETER FOR EACH ACTIVITY 
GHOST THEPACTIVITY CONTAINS AN ALTERNATE EXIT. 
ADO INO RCATES NO NEXT ACTIVITY 
ACTAVETY AL 
ie 
2 
ACGEIVI LY G2 
Be 
3 
AC hI Lele 
is 
m 
PC eV ay wee 
Ej: 
0 
PYPE ONUMBER (OE-TRANSACTION PRIORTTY “CLASSES 
ER 
AL 
TYPE NUMBER OF PARAMETERS FOR EACH TRANSACTION 
ile 
dL, 
LYPE NUMBER OF QUEUES REQUIRED IN MODEL 
fil: 
0) 
FYPE NUMBER OF STORAGES REQUIRED FOLLOWED BY 
CAPACITY SOR TEACE 
bls 
2 5 OO) D'S; (0) 
TYRE NUMBER OCF PACILITIES REQUIRED 
0: 
il 
ARLAANI OT HANCAGCTIONS LO BE PLACED (ON PUTURE 
EVENTS DiST2 AYES OR NO 
NO 
ARE VARIABLES MEAN AND MOD REQUIRED TO GENERATE 
NEWETRANSACTIONS - YES OK NG 
YES 
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TYPE VALUES OF MEAN AND MOD 
O: 


20 6) 
TYRE SENET VALOR OF SIMULATED TIME 
U: 


0 
[TYPE LENGTH OF SIMULATION RUN 
O: 


480 
TYPE MAXIMUM NUMBER OF TRANSACTIONS TO BE TERMINATED 
O: 


1000 
IF ADDITIONAL VARIABLES ARE REQUIRED IN THIS MODEL 
THEY MUST BE SET UP BEFORE RUNNING THE MODEL 


ENTER COMMAND WORD 
RUN 


SIMULATION RUN ENDED AT TIME 476 
NUMBER TRANSACTIONS TERMINATED 24 


ARE STATISTICS REQUIRED? YES OR NO 
YES 

ARE COLUMN TITLES REQUIRED? YES OR NO 
YES 
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MPACI OLE “STATISTICS * 

COL1 FACILITY NUMBER 

COL2 AVERAGE UTILIZATION 
COL3 TOTAL ENTRIES 

COL4 AVERAGE TIME PER ENTRY 
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*STORAGE STATISTICS* 

COL1 STORAGE NUMBER 

COL2 STORAGE CAPACITY 

COL3 CURRENT CONTENTS 

COL4Y MAXIMUM CONTENTS 

COLS AVERAGE CONTENTS 

COL6 AVERAGE UTILIZATION 

COLT TOTAL STORAGE UNITS ENTERED 
COL8 AVERAGE TIME PER STORAGE UNIT 
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Sco here 1.900F1 2.08482 


*QUEUE STATISTICS* 
NO QUEUES IN THIS MODEL 
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END 
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ENTER COMMAND WORD 
NEWMODEL 


TYPE NAMES OF ACTIVITIES TO BE USED IN THIS MODEL 
BACH NAME CONTAINS 6 CHARACTERS SEPARATED BY A BLANK 
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ENDISK PROMES ENDMES 
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TYPE NUMBER OF QUEUES REQUIRED IN MODEL 
O: 


3 
TYPE NUMBER OF STORAGES REQUIRED FOLLOWED BY 
CAPACITY OF EACH 


O: 

al SOO 
TYPE NUMBER OF “FACILITIES REQUIRED 
O: 


a 
ARE ANY TRANSACTIONS TO BE PLACED ON FUTURE 
EVENTS LIST? YES OR NO 
NO 
ARE VARIABLES MEAN AND MOD REQUIRED TO GENERATE 
NEW TRANSACTIONS? YES OR NO 


JARS 
TYPE VALUES OF MEAN AND MOD 
O: 
50 25 
TYPE INIFPLAD* VALUE OF 4SIMULATED?TIME 
O: 


0 
TYPE LENGTH OF SIMULATION RUN 
O: 

3000 
TYPE MAXIMUM NUMBER OF TRANSACTIONS TO BE TERMINATED 
O: 

100 
EPeADDITIONAL-VARFABGES ARE REQUIRED IN THIS MODEL 
THEY MUST BE SET UP BEFORE RUNNING THE MODEL 


ENTER COMMAND WORD 
RUN 


SIMULATION RUN ENDED AT TIME 2995 
NUMBER TRANSACTIONS TERMINATED 47 


ARE STATISTICS REQUIRED? YES OR NO 
YES 

ARE COLUMN TITLES REQUIRED? YES OR NO 
YES 
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*RACIDITY STATISTICS 

COL1 FACILITY NUMBER 

COL2 AVERAGE UTILIZATION 
COL3 TOTAL ENTRIES 

COL4Y AVERAGE TIME PER ENTRY 
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*STORAGE STATISTICS* 

COL1 STORAGE NUMBER 

COL2 STORAGE CAPACITY 

COL3 CURRENT CONTENTS 

COL4 MAXIMUM CONTENTS 

COLS AVERAGE CONTENTS 

COL6 AVERAGE UTILIZATION 

COLT TOTAL STORAGE UNITS ENTERED 
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ENTER COMMAND WORD 
NEWMODEL 


LYPE WAMES OF ACTIVITIES TO BE USED IN THIS MODEL 
EACH NAME CONTAINS 6 CHARACTERS SEPARATED BY A BLANK 


GENRAT ENTER INFORM SPLIT BOOKLI DUPMAC ENDTR2 
DOCLIB FILMQU READER ENDTR1 


HYPE ONE CONFIGURATION PARAMETER FOR BACH ACTIVITY 
POP Paneer py Ly aCONTADNS AN ALTERNATE, EXIT 
ApsOLBI NDE CATE SaNOWNEXE AGTRVITY 
ACPEVETY ¢ at 
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TANCE SD EIEN Pe 
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AGH VILTLY 28 
O: 
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AGRIV ITY 69 
le 
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AGHIVLTL $20 
ile 
li all 
ARG Tele VoIETAY, Sled: 
ae 
6) 
TYPE NUMBER OF TRANSACTION PRIORITY CLASSES 
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TYPE NUMBER OF PARAMETERS FOR FACH TRANSACTION 
Be 

il 
TYPE NUMBER OF QUEUES REQUIRED IN MODEL 
I]: 

5 
TYPE NUMBER OF STORAGES REQUIRED FOLLOWED BY 
CAPACITY OF FACE 


O: 

1 3 
FYPEGNUMBER OF FACEERITIES REQUIRED 
O: 


6 
ARE ANY TRANSACTIONS TO BE PLACED ON FUTURE 
EVENTSTOLat2 YESAQCE TNO 
NO 
ARE VARIABLES MEAN AND MOD REQUIRED TO GENERATE 
NEW TRANSACTIONS? YES OR NO 


YES 
TYPE VALUES OF MEAN AND MOD 
ake 
5 2 
LY PE RINT LIAL VABUE OF. ol MUGAT EO SIME 
Bie 


6) 
TYPE LENGTH OF SIMULATION RUN 
cs 

1000 
TYPE MAXIMUM NUMBER OF TRANSACTIONS TO BE TERMINATED 
UO: 

100 
IF ADDITIONAL VARIABLES ARE REQUIRED IN THIS MODEL 
THEY MUST BE SET UP BEFORE RUNNING THE MODEL 


ENTER COMMAND WORD 
RUN 


SIMULATION RUN ENDED AT TIME 570 
NUMBER TRANSACTIONS TERMINATED 100 


ARE STATISTICS REQUIRED? YES OR NO 
YES 

ARE COLUMN TITLES REQUIRED? YES OR NO 
YES 
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*PACILITY STATISTICS* 

COL1 FACILITY NUMBER 

COL2 AVERAGE UTILIZATION 
COL3 TOTAL ENTRIES 

COL4 AVERAGE TIME PER ENTRY 


a ey TUL S| 56 7.839 
2 OFF 02.6 16 20,10 
3 75189 5 41.2 

4 0.8268 26 ab f 5 cshal 
5 0.4138 iba 20003 
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*STORAGE STATISTICS* 

COL1 STORAGE NUMBER 

COL2 STORAGE CAPACITY 

COL3 CURRENT CONTENTS 

COL4Y MAXIMUM CONTENTS 

COLS AVERAGE CONTENTS 

COL6 AVERAGE UTILIZATION 

COLT TOTAL STORAGE UNITS ENTERED 
COL8 AVERAGE TIME PER STORAGE UNIT 


a 3 2 3 On422.9 
OReL td. i SHEN ea) 


*QUEUE STATISTICS* 

COL1 QUEUE NUMBER 

COL2 CURRENT CONTENTS 

COL3 MAXIMUM CONTENTS 

COL4Y AVERAGE CONTENTS 

COLS TOTAL ENTRIES 

COL6 ZERO DELAY ENTRIES 

COLT AVERAGE TIME PER ENTRY 

COL8 AVERAGE TIME PER ENTRY EXCLUDING 
ZERO DELAY ENTRIES 
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END 


SIMULATION COMPLETED 
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